SYSTEM ONLINE — FYP PROJECT 2024 ACADEMIC USE ONLY
UNCLASSIFIED // ACADEMIC

Small-Scale BALLISTIC MISSILE Development System

[ Final Year Project — Integrated Design, Propulsion & Guidance Research ]

An academic-grade engineering research project covering structural design, solid propellant propulsion, flight dynamics, guidance systems, and trajectory simulation for a small-scale ballistic missile prototype.

VIEW MISSION STAGES PROJECT BRIEF
0 Stages
|
0 Est. Days
|
0 Engineers
|
0 Tools
SCROLL
PROJECT OVERVIEW

Small-Scale Ballistic Missile
Final Year Project

A comprehensive academic research project simulating real-world ballistic missile development at a reduced academic scale — covering every engineering discipline from propellant chemistry to guidance algorithms.

Propulsion System

Solid-propellant rocket motor design, thrust curve analysis, and nozzle geometry optimization for controlled burn performance.

Structural Design

Lightweight airframe engineering using composite materials, aeroelastic analysis, and load-bearing optimization for flight loads.

Guidance & Control

IMU-based inertial guidance, fin-actuated control surfaces, and PID flight controller for trajectory correction mid-flight.

Trajectory Analysis

6-DOF flight simulation, ballistic flight path prediction, range/altitude trade studies, and impact point estimation.

CFD / Aerodynamics

ANSYS-based airflow simulation to quantify drag, stability derivatives, and shock wave behavior at transonic speeds.

Avionics & Sensors

On-board data acquisition, GPS/INS integration, telemetry downlink, and flight computer hardware selection and integration.

ENGINEERING TOOLS:
SolidWorks ANSYS Fluent MATLAB / Simulink OpenRocket Python / SciPy Fusion 360 Arduino / STM32 RASAero II
MISSION FLOWCHART

Development Stage Sequence

Click any stage to expand classified technical details for that phase.

Design
CONCEPT DESIGN PHASE

Conceptual Design & Mission Requirements

Every ballistic missile project begins by defining the mission envelope — target range, payload mass, apogee altitude, and acceptable impact dispersion. These requirements drive all downstream design decisions.

▶
Mission Definition: Target range (100–500 m for small-scale), maximum altitude (150–500 m), payload type (data logger or dummy mass), and recovery system requirements.
▶
Geometric Sizing: Fineness ratio (L/D = 10–15) selected to balance drag and structural length. Nose cone profile chosen from Von Kármán, ogive, or conical geometry based on Mach target.
▶
Mass Budget: Component-level mass estimation allocating weight to airframe, motor, avionics, payload, and recovery. Target launch mass 0.5–3 kg for small-scale prototype.
▶
Stability Analysis: Static stability margin (Cp − Cg ≥ 1 caliber) verified using Barrowman equations before proceeding to detailed design.
▶
Safety & Legal: PCAA/SUPARCO requirements for amateur rocketry in Pakistan, range safety plan, and academic institution (HEC-recognized university) approval process documented at this stage.
Barrowman EqnsMass BudgetOpenRocketPCAA Compliance
Propulsion
PROPULSION PHASE

Solid Propellant Propulsion System

The propulsion system is the heart of any ballistic missile. For small-scale academic projects, commercial off-the-shelf (COTS) solid rocket motors (APCP) are used, or custom grain geometries are studied analytically.

▶
Motor Selection: Aerotech/Cesaroni APCP motors (class F–J) selected based on required total impulse. Thrust-to-weight ratio target ≥ 5:1 at launch.
▶
Nozzle Design: Convergent-divergent (De Laval) nozzle geometry designed using isentropic flow equations. Throat diameter and exit area ratio optimized for sea-level expansion.
▶
Grain Geometry: BATES (Ballistic Test and Evaluation System) grain configuration studied. Burn surface area vs. regression rate analysis produces predicted thrust curve.
▶
Specific Impulse: Isp target of 180–220 s for APCP compositions. Chemical equilibrium calculations using NASA CEA code validate propellant performance.
▶
Motor Casing: Aluminum or filament-wound CFRP casing sized for chamber pressure (typically 30–80 bar) with safety factor ≥ 4.
APCP MotorsDe Laval NozzleNASA CEABATES GrainIsp
CAD Model
CAD MODELLING PHASE

3D Structural & CAD Modelling

Full parametric 3D models are constructed in SolidWorks or Fusion 360, capturing every structural component of the missile from nose tip to nozzle exit plane.

▶
Airframe Tube: Main body tube modelled as a thin-walled cylinder with internal bulkheads, centering rings, and motor retention components all parametrically defined.
▶
Nose Cone: Ogive or Von Kármán profile generated using parametric equations and 3D-printed in PETG or machined from aluminum depending on Mach number requirements.
▶
Fin Assembly: Trapezoidal or clipped-delta fins with 3 or 4-fin configuration. Root chord, tip chord, and sweep angle optimized using RASAero II stability output.
▶
Recovery Bay: Parachute deployment section modelled with ejection charge housing, shock cord attachment points, and nose cone separation mechanism.
▶
FEA Integration: SolidWorks Simulation performs stress analysis under max-Q dynamic pressure load, motor thrust load, and landing impact scenario.
SolidWorksFusion 360FEA3D PrintingRASAero II
CFD
CFD ANALYSIS PHASE

Computational Fluid Dynamics Analysis

ANSYS Fluent is used to simulate compressible airflow around the missile body from launch through max-Q, generating aerodynamic coefficients used in the 6-DOF trajectory simulation.

▶
Mach Sweep: Steady-state simulations at Mach 0.2 to 1.5 generate Cd, CL, and Cm as functions of angle-of-attack (±15°) and Mach number.
▶
Shock Analysis: Transonic and supersonic cases capture bow shock formation, shock-boundary layer interaction at fins, and base pressure drag contribution.
▶
Turbulence Model: k-ω SST two-equation model employed for accurate boundary layer transition prediction on fin leading edges and nose cone.
▶
Thermal Analysis: Coupled thermal-structural solver quantifies aerodynamic heating at stagnation point and leading edges using recovery temperature and heat transfer coefficients.
▶
Grid Independence: Mesh refinement study conducted with y+ ≤ 1 near walls to ensure simulation accuracy independent of grid density.
ANSYS Fluentk-ω SSTMach SweepMax-QShock Analysis
GNC
GNC SYSTEMS PHASE

Guidance, Navigation & Control (GNC)

The GNC system is the brain of the missile, processing sensor data in real-time to maintain the desired trajectory and orientation throughout the ballistic flight phase.

▶
IMU Selection: MPU-6050 or BNO055 9-DOF inertial measurement unit provides 3-axis accelerometer, gyroscope, and magnetometer data at 100–200 Hz sample rate.
▶
Flight Computer: STM32 or Teensy 4.1 microcontroller runs real-time GNC algorithms. Barometric altimeter (BMP388) and GPS module provide supplementary state estimation.
▶
PID Controller: Three-axis PID attitude controller drives four servo-actuated fins (or thrust vector control carbon vanes) to correct pitch, yaw, and roll deviations.
▶
Kalman Filter: Extended Kalman Filter (EKF) fuses IMU, barometer, and GPS data for accurate state estimation (position, velocity, attitude) under noisy sensor conditions.
▶
Telemetry: 433 MHz LoRa or 900 MHz SiK radio transmits live flight data (altitude, velocity, attitude) to ground station for real-time monitoring and post-flight analysis.
STM32PID ControlKalman FilterIMULoRa Telemetry
Trajectory
TRAJECTORY PHASE

Trajectory Simulation & Flight Analysis

High-fidelity trajectory simulation integrates the full aerodynamic database, propulsion thrust curve, and mass properties to predict the complete flight from ignition to impact.

▶
6-DOF Model: Full six-degree-of-freedom equations of motion coded in MATLAB/Simulink, accounting for translational and rotational dynamics simultaneously.
▶
OpenRocket Validation: OpenRocket provides rapid first-order trajectory prediction, used to cross-validate the higher-fidelity MATLAB 6-DOF model outputs.
▶
Launch Angle Study: Parametric sweep of launch elevation angle (45°–85°) vs. range and altitude output generates optimal launch angle for target engagement.
▶
Wind Effects: Monte Carlo analysis with randomized wind speed and direction inputs produces statistical impact dispersion ellipses (CEP — Circular Error Probable).
▶
Recovery Simulation: Parachute deployment altitude (at apogee), descent rate, and drift distance under wind modelled to ensure safe landing zone prediction.
MATLABSimulink6-DOFOpenRocketMonte CarloCEP
Materials
MATERIALS & MFG PHASE

Materials Selection & Manufacturing

Material choice directly determines structural performance, weight, and manufacturability. The following materials are evaluated for each subsystem of the small-scale missile.

AIRFRAME TUBE Carbon Fiber / Fiberglass (G12)

High-modulus CFRP or phenolic tubes. CFRP: 1.6 g/cm³, 70 GPa modulus. Fiberglass G12 more cost-effective for subsonic flights.

NOSE CONE 3D-Printed PETG / Machined Aluminum

PETG 3D-printed nose cones sufficient for M<0.8. Aluminum 6061-T6 nose machined for transonic/supersonic flights requiring thermal resistance.

FIN SET Plywood G10 / Aluminum 6061

Aircraft plywood (6 mm) for low-cost subsonic designs. G10 fiberglass or 1.5 mm aluminum 6061 sheet for structural fins at higher velocities.

MOTOR CASING Aluminum 7075-T6 / CFRP Overwrap

Aluminum 7075-T6 for commercial motor casings. Filament-wound CFRP provides 40% weight saving for custom motor design at higher chamber pressures.

CFRPAluminum 6061G10/G12PETG7075-T6
Testing
TESTING & LAUNCH PHASE

Testing, Launch Operations & Recovery

Ground testing and live launch are the culmination of all prior engineering work. A strict safety and test protocol is followed to protect personnel, equipment, and comply with applicable regulations.

▶
Static Motor Test: Motor fired on a thrust stand instrumented with a load cell and pressure transducer. Measured thrust curve compared against manufacturer data and simulation.
▶
Avionics Bench Test: Full GNC system powered on bench with simulated IMU data injected to verify controller response, data logging, and parachute ejection logic.
▶
Pre-Launch Checklist: 25-point launch readiness review covering structural, propulsion, avionics, safety, range clearance, and weather window checks.
▶
Live Launch: Conducted at a designated open field with university/SUPARCO approval in Pakistan. Launch rail angle set per trajectory simulation optimal angle. Telemetry monitored from ground station.
▶
Post-Flight Analysis: Logged flight data (altitude, velocity, attitude) compared against simulation predictions. Drag coefficient, peak altitude, and apogee time used for model validation and iteration.
Static FireThrust StandSUPARCO/UniversityTelemetryPost-Flight Analysis
TECHNICAL REFERENCE

Ballistic Missile Classification & Types

Understanding the full spectrum of ballistic missile systems — from tactical short-range to strategic intercontinental platforms.

Ballistic missiles are categorized primarily by their operational range and launch platform. This classification determines the propulsion requirements, guidance complexity, and strategic role of each missile system.

Range-Based Classification

Type Abbreviation Range (km) Range (miles) Examples
Tactical Ballistic Missile TBM < 300 km < 190 mi ATACMS, Tochka, Iskander-M
Short-Range Ballistic Missile SRBM 300 – 1,000 km 190 – 620 mi Scud-B, Prithvi-I, Fateh-110
Medium-Range Ballistic Missile MRBM 1,000 – 3,500 km 620 – 2,170 mi Agni-II, Shaheen-II, DF-21
Intermediate-Range Ballistic Missile IRBM 3,500 – 5,500 km 2,200 – 3,400 mi Agni-III, Shaheen-III, DF-26
Intercontinental Ballistic Missile ICBM > 5,500 km > 3,400 mi Minuteman-III, RS-24 Yars, DF-41
Submarine-Launched Ballistic Missile SLBM Variable Variable Trident-II, K-4, R-30 Bulava

Source: Defense Intelligence Ballistic Missile Analysis Committee (2017), Praeger Security International

Boost Phase

Powered flight duration from ignition to fuel exhaustion. Multi-stage rocket propulsion with internal computers maintaining preprogrammed trajectory. Stage separation occurs during this phase.

  • Duration: Tenths of seconds to several minutes
  • Ends when all fuel is exhausted
  • Multi-stage separation sequence
  • Maximum acceleration and dynamic pressure

Mid-Course Phase

Longest trajectory phase with missile in free flight (no thrust). May travel through space at extreme altitudes before re-entry begins.

  • Apogee: ~4,500 km for ICBMs
  • Velocity: 7.5–10 km/s (4–5 nautical miles/s)
  • Consists of warhead/payload with possible countermeasures
  • Potential deployment of MIRV and decoys

Terminal/Re-entry Phase

Begins at atmospheric re-entry where aerodynamic forces become significant. Warhead experiences extreme heating and deceleration.

  • Re-entry velocity: 6–8 km/s (22,000–29,000 km/h)
  • Extreme aerodynamic heating (up to 2,000°C)
  • Atmospheric drag significantly affects trajectory
  • Final impact point determination

Throw-Weight

Measure of effective payload weight (kg or tonnes) including warheads, re-entry vehicles, dispensing mechanisms, and penetration aids. Excludes booster rockets and fuel. Used as classification criterion in arms control treaties (SALT).

Heaviest: Russian SS-18, Chinese CSS-4

MIRV Technology

Multiple Independently Targetable Re-entry Vehicles allow a single missile to carry multiple warheads, each capable of striking different targets. Post-boost vehicle (PBV) or "bus" maneuvers to release RVs + decoys + chaff.

First MIRV ICBM: Minuteman III (1970)

Trajectory Variants

Different flight profiles serve different strategic purposes. Minimum-energy maximizes payload, depressed trajectory reduces flight time, and lofted trajectory is used for testing.

Depressed: Evades ABM systems

Hypersonic Flight

Modern ballistic missiles reach Mach 5+ during re-entry. Hypersonic glide vehicles (HGV) perform atmospheric maneuvers with real-time aerodynamic calculations, distinguishing from traditional ballistic trajectory.

Examples: Kinzhal, Oreshnik, Avangard
ENGINEERING GUIDE

Missile Structural Design & Manufacturing

Complete technical guide to designing and building a small-scale ballistic missile — from airframe geometry to final assembly.

1. Airframe / Body Tube Design

Geometry & Sizing

  • Diameter: 75–150 mm for small-scale (based on motor size)
  • Length: Determined by fineness ratio (L/D = 10–15)
  • Wall Thickness: 2–4 mm for composite tubes, 1.5–3 mm for aluminum
  • Material: CFRP (carbon fiber), fiberglass G12, or aluminum 6061-T6
  • Surface Finish: Smooth exterior to minimize skin friction drag

Structural Considerations

  • Max-Q Resistance: Designed to withstand maximum dynamic pressure without buckling
  • Motor Thrust Load: Bulkheads transfer thrust to airframe
  • Landing Impact: Reinforced sections at recovery attachment points
  • Thermal Protection: Internal insulation near motor section
Design Formula: Fineness Ratio = Length / Diameter (Optimal: 10–15)
Critical Buckling Pressure: P_cr = (2×E×t³) / (D³×(1-ν²))

2. Nose Cone Design

Profile Type Equation Best For Drag Characteristics
Conical y = (R/L)×x Subsonic (M < 0.8) Higher drag, simple to manufacture
Tangent Ogive y = R - √(R² - x²) Transonic (M = 0.8–1.2) Good balance, popular choice
Von Kármán y = f(x, M, L/D) Supersonic (M > 1.2) Minimum wave drag, optimal
Parabolic y = (R/L²)×x² Low supersonic Simple, moderate drag

Manufacturing Options

  • 3D Printing: PETG/PLA for M < 0.8, infill ≥ 40%
  • CNC Machining: Aluminum 6061-T6 for transonic/supersonic
  • Composite Layup: Carbon fiber over foam core
  • Wood Turning: Hardwood for low-speed prototypes

Attachment Methods

  • Shoulder Fit: Press-fit with friction retention
  • Screw Attachment: Internal threads with retention screw
  • Bayonet Mount: Twist-lock for quick release
  • Shear Pin: For parachute deployment separation

3. Fin Design & Attachment

Fin Geometry Parameters

  • Number of Fins: 3-fin (lower drag) or 4-fin (better stability)
  • Planform: Trapezoidal, clipped-delta, or elliptical
  • Root Chord: 25–40% of body length
  • Tip Chord: 40–60% of root chord
  • Semi-span: 15–25% of body diameter
  • Sweep Angle: 30–45° for transonic/supersonic
  • Thickness: 3–6 mm (airfoil or flat plate)

Materials & Construction

  • Aircraft Plywood: 3–6 mm birch, epoxy-coated (subsonic)
  • G10 Fiberglass: 2–4 mm, excellent strength-to-weight
  • Aluminum 6061: 1.5–3 mm sheet, machined or waterjet
  • Carbon Fiber: Pre-preg layup, autoclave cured (high performance)
  • 3D Printed: PETG/ABS with ≥60% infill (prototypes)
Critical Design Check:

Ensure static stability margin (Cp - Cg) ≥ 1 caliber (body diameter). Use Barrowman equations or RASAero II for calculation. Fin flutter velocity must exceed maximum expected speed by 20% safety margin.

Fin Attachment Methods

Through-Wall Mounting:

Fin root passes through airframe wall, bonded with epoxy. Strongest method, requires precise cutting.

Surface Mount:

Fin bonded to exterior surface. Simpler but weaker, suitable for small fins.

Fin Can:

Fins attached to separate tube that slides into main airframe. Modular design.

Integrated:

3D printed or molded as single piece with airframe. No bonding required.

4. Motor Mount Assembly

Motor Tube Specifications

  • Material: Aluminum 6061-T6 or phenolic tube
  • Wall Thickness: 3–5 mm (depends on chamber pressure)
  • Length: Motor length + 20–30 mm clearance
  • Inner Diameter: Motor outer diameter + 0.5–1 mm slip fit
  • Thermal Protection: Ceramic coating or insulating liner

Centering Rings

  • Material: 6–10 mm plywood, G10, or aluminum
  • Quantity: 2 rings (forward and aft)
  • Function: Center motor tube in airframe, transfer thrust
  • Bonding: Epoxy to both motor tube and airframe interior
  • Vent Holes: 3–4 holes for pressure equalization

Thrust Transfer

  • Aft Ring: Primary thrust bearing surface
  • Retainer: Screw or clip prevents motor ejection
  • Forward Ring: Locates motor, minimal thrust load
  • Engine Hook: For commercial motor retention

Nozzle Clearance

  • Exit Distance: Nozzle extends 10–20 mm aft of airframe
  • Heat Shield: Stainless steel or ceramic at aft end
  • Exhaust Path: Unobstructed flow, no impingement

5. Bulkheads & Recovery System

Bulkhead Design

  • Forward Bulkhead: Separates payload/avionics from recovery bay
  • Aft Bulkhead: Separates recovery bay from motor section
  • Material: 6–12 mm plywood or aluminum plate
  • Attachment: Epoxy bonded + screws for redundancy
  • Wire Pass-Through: Grommet-protected holes for avionics wiring

Recovery System Components

  • Parachute: Drogue (small, high-altitude) + main (large, landing)
  • Shock Cord: Kevlar or nylon webbing, 2–3× rocket length
  • Ejection Charge: Black powder charge (1–3 grams) for separation
  • Piston/Plate: Pushes nose cone out without damaging parachute
  • Deployment Bag: Protects parachute during ejection
Recovery Sequence:
  1. Apogee detected by flight computer (accelerometer + barometer)
  2. Drogue ejection charge fires → nose cone separates → drogue deploys
  3. Descent to ~300 m altitude
  4. Main ejection charge fires → drogue jettisoned → main parachute deploys
  5. Controlled descent at 5–7 m/s to landing zone

6. Step-by-Step Assembly Procedure

Phase 1: Airframe Preparation

  1. Cut body tube to final length using fine-tooth saw or tube cutter
  2. Sand cut edges smooth, deburr interior and exterior
  3. Mark fin slot locations using alignment jig (equal spacing: 120° for 3-fin, 90° for 4-fin)
  4. Cut fin slots using router or Dremel with cutting guide (depth = 1.5× fin root thickness)
  5. Drill vent holes near forward and aft ends (4× 6 mm holes for pressure equalization)
  6. Clean interior surface with isopropyl alcohol for bonding preparation

Phase 2: Motor Mount Assembly

  1. Cut motor tube to length (motor + 25 mm)
  2. Cut centering rings from plywood/G10 (outer diameter = airframe ID, inner = motor tube OD)
  3. Test fit motor tube through centering rings (should slide with slight friction)
  4. Apply epoxy to motor tube, slide on centering rings, wipe excess
  5. Allow epoxy to cure 24 hours at room temperature
  6. Install motor retainer (screw threads or engine hook) at aft end
  7. Slide motor mount assembly into airframe, bond aft centering ring to airframe

Phase 3: Fin Installation

  1. Cut fins from material using template (ensure identical geometry)
  2. Sand fins to airfoil shape (rounded leading edge, sharp trailing edge)
  3. Dry-fit fins into slots (should fit snugly without gaps)
  4. Apply epoxy to fin root and slot walls
  5. Insert fins, ensure perpendicular alignment (use square tool)
  6. Apply fillet bead along fin-airframe joint (epoxy + microballoons)
  7. Allow 24 hours cure time, sand fillets smooth
  8. Prime and paint fins if desired (multiple light coats)

Phase 4: Bulkheads & Recovery

  1. Cut bulkheads from plywood/aluminum (diameter = airframe ID - 0.5 mm)
  2. Drill wire pass-through holes, install grommets
  3. Bond aft bulkhead above motor mount (epoxy + screws)
  4. Install shock cord anchor point on forward bulkhead
  5. Assemble parachute, shock cord, and deployment bag
  6. Load black powder ejection charge in deployment well
  7. Thread avionics wiring through bulkhead pass-throughs

Phase 5: Nose Cone & Final Assembly

  1. Machine or 3D print nose cone to specification
  2. Sand nose cone shoulder for smooth fit (no binding, no wobble)
  3. Install retention screw or bayonet mount
  4. Attach shock cord to nose cone base
  5. Install avionics bay (flight computer, battery, sensors)
  6. Pack recovery system (deployment bag → parachute → shock cord)
  7. Final balance check: locate CG with motor installed
  8. Apply final paint/decal scheme (lightweight paint only)
AUTOMATION SYSTEMS

Missile Automation & Flight Control

Complete guide to automating a ballistic missile — flight computer architecture, sensor fusion, PID control, and real-time guidance algorithms.

1. Flight Computer Hardware Architecture

Microcontroller Options

STM32F4/F7 168–216 MHz, FPU, DSP, 512KB–2MB Flash
Teensy 4.1 600 MHz, 1MB Flash, excellent for rapid prototyping
Pixhawk (Cube) Dual IMU, redundant sensors, aerospace-grade
Arduino Due 84 MHz ARM Cortex-M3, beginner-friendly

Key Requirements

  • Clock Speed: ≥ 100 MHz for real-time control loops
  • FPU: Hardware floating-point for fast math operations
  • ADC: 12-bit+ for sensor readings
  • Timers: Multiple hardware timers for PWM generation
  • Communication: UART, SPI, I2C for sensor interfacing
  • Memory: ≥ 256 KB Flash, ≥ 64 KB RAM

Recommended: STM32F405/F446

CPU:

ARM Cortex-M4 @ 168 MHz with FPU and DSP instructions

Memory:

192 KB SRAM, 1 MB Flash, external SPI Flash support

Peripherals:

3× SPI, 3× I2C, 6× UART, 12× timers, 3× ADC (12-bit)

Why STM32?

Industry standard, excellent HAL library, low power, widely documented

2. Sensor Suite & Interfacing

IMU (Inertial Measurement Unit)

Primary attitude and acceleration sensing

MPU-60506-DOF (gyro+accel), I2C, budget option
BNO0559-DOF with sensor fusion, UART/I2C
ICM-209489-DOF, low noise, SPI/I2C
BMI088High-g accel, aerospace-grade

Interface: I2C (400 kHz) or SPI (10 MHz)

Barometric Altimeter

Altitude and vertical velocity measurement

BMP388±0.5 m accuracy, I2C/SPI
MS5611±10 cm accuracy, high-speed
LPS25HLow power, ±0.3 m accuracy

Interface: I2C or SPI

GPS Module (Optional)

Position and velocity for long-range flights

UBlox NEO-M8N10 Hz, 2.5 m CEP
BN-880Built-in compass, 10 Hz
F9P RTKCentimeter-level accuracy

Interface: UART (9600–115200 baud)

Additional Sensors

  • Magnetometer: HMC5883L, QMC5883L (heading reference)
  • Current Sensor: INA219 (battery monitoring)
  • Voltage Divider: Battery voltage monitoring
  • Temperature: Internal IMU or external thermistor

Sensor Fusion Architecture

IMU (Gyro+Accel)

100–200 Hz raw data

Barometer

50–100 Hz altitude

GPS (Optional)

10 Hz position

Kalman Filter

Fused state estimate

3. PID Flight Control System

PID Controller Fundamentals

PID (Proportional-Integral-Derivative) control is the industry standard for missile attitude control. It calculates the error between desired and actual attitude, then applies corrective action.

output = Kp×error
+ Ki×∫error×dt
+ Kd×d(error)/dt
  • Kp (Proportional): Response to current error
  • Ki (Integral): Eliminates steady-state error
  • Kd (Derivative): Dampens oscillations, predicts future

Three-Axis Control

Missile attitude is controlled along three rotational axes. Each axis has an independent PID controller.

Pitch Nose up/down rotation
Yaw Nose left/right rotation
Roll Clockwise/counterclockwise rotation
Control Loop Rate: 500–1000 Hz for stable flight

Fin Actuation (Servo Motors)

  • Servo Type: Digital high-speed servos (0.06–0.08 sec/60°)
  • Torque: ≥ 5 kg·cm for small fins, ≥ 10 kg·cm for larger
  • Examples: BMS-2275, MKS DS95i, Hitec D85MG
  • Control Signal: PWM (50 Hz, 1000–2000 μs pulse width)
  • Power: 5–6 V BEC (Battery Eliminator Circuit)
  • Mounting: Direct horn attachment or pushrod linkage

Thrust Vector Control (TVC)

  • Method: Gimbaled nozzle with 2–4 actuators
  • Actuators: High-torque servos or linear actuators
  • Deflection: ±3–5° vectoring angle
  • Response: Must be faster than aerodynamic control
  • Advantage: Control at low speeds (before fins are effective)
  • Complexity: Higher than fin control, requires robust mechanism

PID Tuning Procedure

Step 1: Initialize

Set Ki = 0, Kd = 0. Start with low Kp (0.1–0.5).

Step 2: Tune Kp

Increase Kp until system oscillates, then reduce by 30%.

Step 3: Tune Kd

Add derivative gain to dampen oscillations. Watch for noise amplification.

Step 4: Tune Ki

Add integral gain to eliminate steady-state error. Use sparingly.

Step 5: Flight Test

Fine-tune gains based on actual flight performance. Log data for analysis.

Typical Values

Kp: 0.5–2.0, Ki: 0.001–0.01, Kd: 0.05–0.2 (axis-dependent)

4. Guidance & Navigation Algorithms

Extended Kalman Filter (EKF)

The EKF fuses IMU, barometer, and GPS data to estimate the missile's state (position, velocity, attitude) with minimal error.

Predict Step:

x̂ₖ⁻ = F×x̂ₖ₋₁ + B×uₖ
Pₖ⁻ = F×Pₖ₋₁×Fᵀ + Q

Update Step:

Kₖ = Pₖ⁻×Hᵀ×(H×Pₖ⁻×Hᵀ + R)⁻¹
x̂ₖ = x̂ₖ⁻ + Kₖ×(zₖ - H×x̂ₖ⁻)

Flight State Machine

The flight computer operates as a finite state machine, transitioning through predefined flight phases.

  1. IDLE: Pre-launch, system checks
  2. LAUNCH DETECT: Acceleration threshold trigger
  3. BOOST: Powered ascent, active control
  4. COAST: Ballistic ascent to apogee
  5. APOGEE: Drogue parachute deployment
  6. DESCENT: Main parachute deployment
  7. LANDING: Touchdown detected, data save

Launch Detection Algorithm

Launch is detected by monitoring accelerometer data for a sudden increase in acceleration exceeding the gravity threshold.

// Launch detection threshold
const float LAUNCH_ACCEL = 2.5; // g

if (accel_magnitude > LAUNCH_ACCEL) {
  flight_state = BOOST;
  launch_time = millis();
}
Apogee Detection

Apogee is detected when vertical velocity crosses zero (positive to negative) or barometric altitude reaches maximum.

if (vertical_velocity < -0.5 &&
  prev_velocity >= 0) {
  deploy_drogue();
  flight_state = APOGEE;
}

5. Software Architecture & Code Structure

Main Loop Structure

void loop() {
  // Read sensors (1000 Hz)
  readIMU();
  readBarometer();
  readGPS();

  // Sensor fusion (500 Hz)
  kalmanUpdate();

  // Flight state machine
  updateFlightState();

  // PID control (500-1000 Hz)
  pidCompute();
  updateServos();

  // Data logging (50 Hz)
  logFlightData();
  telemetryTransmit();
}

Key Software Modules

  • sensor_drivers.cpp: IMU, barometer, GPS drivers
  • kalman_filter.cpp: EKF implementation
  • pid_controller.cpp: 3-axis PID control
  • flight_state.cpp: State machine logic
  • servo_output.cpp: PWM generation
  • data_logger.cpp: SD card logging
  • telemetry.cpp: LoRa radio communication

Development Environment

IDE:

STM32CubeIDE, PlatformIO, or Arduino IDE

Debugging:

ST-Link V2 debugger, serial printf, logic analyzer

Version Control:

Git + GitHub for code management

Simulation:

MATLAB/Simulink for control design

6. Telemetry & Ground Station

Radio Communication

LoRa (433/915 MHz) Long range (5–10 km), low data rate
SiK Radio (915 MHz) 300 m–2 km, MAVLink protocol
nRF24L01 (2.4 GHz) Short range (100 m), high data rate
XBee (2.4 GHz) Medium range (1–2 km), reliable

Telemetry Data

  • Altitude (barometric + GPS)
  • Velocity (vertical + horizontal)
  • Attitude (pitch, roll, yaw)
  • Acceleration (3-axis)
  • GPS coordinates
  • Battery voltage
  • Flight state
  • Timestamp

Ground Station Software Options

Mission Planner

Windows-based, MAVLink protocol, full-featured

QGroundControl

Cross-platform, modern UI, MAVLink

Custom Python GUI

PyQt/Matplotlib, fully customizable

Arduino Serial Plotter

Simple debugging, basic visualization

PROJECT TIMELINE

Mission Schedule

Estimated duration for each development phase. Total project span: approximately 62 working days.

01
Conceptual Design 8 Days

Mission spec, geometry sizing, stability checks, safety approvals

02
Propulsion Design 10 Days

Motor selection, nozzle sizing, grain analysis, Isp calculations

03
3D Modelling 8 Days

Full CAD assembly, FEA stress analysis, manufacturing drawings

04
CFD Simulation 12 Days

ANSYS setup, Mach sweep runs, post-processing, aero database

05
GNC Development 12 Days

Flight computer coding, PID tuning, Kalman filter, bench testing

06
Manufacturing 8 Days

Airframe assembly, fin bonding, avionics integration, recovery rigging

07
Testing & Launch 4 Days

Static fire, pre-launch checklist, live launch, data recovery

Total Estimated Duration: 62 Days
BUDGET ESTIMATION

Academic Project Cost Breakdown

Estimated procurement and operational costs for a small-scale ballistic missile FYP prototype in Pakistan. SolidWorks & ANSYS costs excluded (already owned).

# Item Purpose Est. Cost (PKR)
01 Rocket Motors (APCP)Aerotech / Cesaroni F–J class motors (3–5 units), imported via Daraz/AliExpressPKR 35,000 – 85,000
02 Airframe MaterialsCFRP tubes, fiberglass body, centering rings (local composites market, Lahore/Karachi)PKR 22,000 – 56,000
03 3D PrintingNose cone, fin can, payload bay components (PETG/PLA) — university lab or local 3D print shopPKR 11,000 – 28,000
04 Flight ComputerSTM32 / Teensy 4.1, PCB fabrication (JLCPCB/local), connectorsPKR 17,000 – 42,000
05 Sensors & AvionicsIMU, barometer, GPS, LoRa radio, servos (Daraz/AliExpress)PKR 22,000 – 50,000
06 Recovery SystemDrogue + main parachute, shock cord, ejection charge (custom-sewn locally)PKR 11,000 – 28,000
07 CAD/Simulation SoftwareSolidWorks & ANSYS — Already owned (₨0); MATLAB student lic. + free toolsPKR 0 – 8,000
08 Test EquipmentThrust stand, load cell, data logger, pressure gauge (Hall Road / online)PKR 28,000 – 70,000
09 Safety & Launch FeesUniversity ground booking, safety equipment, local transportPKR 22,000 – 56,000
10 Documentation & PrintingReport printing, poster, spiral binding, HEC submissionPKR 8,000 – 22,000
TOTAL ESTIMATED BUDGET PKR 176,000 – 445,000
MISSION TEAM

Required Engineering Personnel

Recommended 5-person team structure for full-scope project execution.

ENG-01

Structural / Airframe Engineer

Leads all CAD modelling, FEA analysis, material selection, manufacturing coordination, and prototype assembly.

SolidWorksFEACFRP
ENG-02

Propulsion Engineer

Manages motor selection, nozzle design, static fire testing, thrust analysis, and propulsion system integration.

NASA CEAAPCPNozzle Design
ENG-03

CFD / Aerodynamics Engineer

Operates ANSYS Fluent for all simulation campaigns, post-processes aero data, and validates against analytical predictions.

ANSYSCFDMeshing
ENG-04

GNC / Avionics Engineer

Develops flight computer firmware, GNC algorithms, sensor integration, telemetry system, and conducts bench testing.

STM32Kalman FilterPID
ENG-05

Flight Dynamics / Trajectory Engineer

Builds and validates 6-DOF MATLAB model, conducts launch angle optimization, Monte Carlo analysis, and post-flight data correlation.

MATLAB6-DOFOpenRocket
ACADEMIC ASSESSMENT

FYP Suitability Analysis

HIGHLY SUITABLE

Recommended for Final Year Project

This project comprehensively satisfies the academic requirements of a mechanical, aerospace, or mechatronics engineering FYP at any HEC-recognized Pakistani university. It integrates six distinct engineering disciplines, produces measurable physical deliverables, and has genuine research novelty through its AI-assisted GNC development.

90%

FYP Suitability Score

Multi-Disciplinary Integration

Combines propulsion, structures, aerodynamics, GNC, and software — meeting broad assessment criteria across multiple engineering modules simultaneously.

Industry Skill Development

Hands-on ANSYS, SolidWorks, STM32 programming, and Python/MATLAB flight simulation directly maps to Pakistan's defence (NESCOM, AWC, NDC), space (SUPARCO), and aerospace industry requirements.

Concrete Physical Deliverable

A working prototype (even static display) with full simulation suite, technical report, and live telemetry data provides compelling assessment evidence.

Research & Publication Potential

ML-based aerodynamic surrogate modeling and low-cost GNC system design are active research areas with potential for student journal submissions.

Fully Compliant & Safe

Commercial APCP motors and supervised launches with university/SUPARCO approval ensure the project remains within academic safety standards and Pakistani regulations.

Measurable Performance Metrics

Apogee altitude accuracy (sim vs. flight), CEP, drag coefficient correlation, and GNC stability margins provide quantitative grading benchmarks.

ACADEMIC REFERENCES

Research Sources & Bibliography

Comprehensive list of technical references, academic papers, and authoritative sources used for this ballistic missile research project.

Ballistic Missiles

  • [1] Defense Intelligence Ballistic Missile Analysis Committee. Ballistic Missiles of the World. Praeger Security International, 2017.
  • [2] Futter, Andrew. Ballistic Missile Defence and US National Security Policy. Routledge, 2013.
  • [3] Neufeld, Jacob. The Development of Ballistic Missiles in the United States Air Force, 1945-1960. Smithsonian Institution Press, 1990.
  • [4] Wikipedia Contributors. "Ballistic Missile." Wikipedia, The Free Encyclopedia. Wikimedia Foundation.

Rocket Propulsion

  • [5] NASA Chemical Equilibrium with Applications (CEA) Code. NASA Glenn Research Center.
  • [6] Wikipedia Contributors. "Solid-propellant rocket." Wikipedia, The Free Encyclopedia. Wikimedia Foundation.
  • [7] Aerotech Inc. "APCP Motor Technology." Technical Documentation, 2023.
  • [8] Cesaroni Technology. "Solid Rocket Motor Design Handbook." 2022.
  • [9] Space Shuttle Solid Rocket Booster Technical Specifications. NASA Historical Archive.

Guidance & Navigation

  • [10] Wikipedia Contributors. "Inertial navigation system." Wikipedia, The Free Encyclopedia. Wikimedia Foundation.
  • [11] Titterton, David H., and John L. Weston. Strapdown Inertial Navigation Technology. IET, 2004.
  • [12] Brown, Robert Grover, and Patrick Y.C. Hwang. Introduction to Random Signals and Applied Kalman Filtering. Wiley, 2012.
  • [13] DARPA Micro-PNT Program. "Timing & Inertial Measurement Unit (TIMU) Technology." 2012.
  • [14] U.S. Army Research Laboratory. "MEMS Sensor Array for Projectile Navigation." 2012.

Aerodynamics & CFD

  • [15] ANSYS Inc. "ANSYS Fluent Theory Guide." 2023.
  • [16] Barrowman, James S. "The Practical Calculation of the Aerodynamic Characteristics of Slender Finned Vehicles." NASA CR-190, 1967.
  • [17] RASAero II. "Rocket Aerodynamic Analysis and Flight Simulation Software." Technical Documentation.
  • [18] Wikipedia Contributors. "Hypersonic flight." Wikipedia, The Free Encyclopedia. Wikimedia Foundation.
  • [19] Sandia National Laboratories. "Hypersonic Materials Testing Program." 2022.

Flight Dynamics

  • [20] OpenRocket Project. "OpenRocket Technical Documentation." 2023.
  • [21] MathWorks. "MATLAB/Simulink 6-DOF Flight Dynamics Modeling." Documentation, 2023.
  • [22] Stengel, Robert F. Flight Dynamics. Princeton University Press, 2004.
  • [23] Zipfel, Peter H. Modeling and Simulation of Aerospace Vehicle Dynamics. AIAA Education Series, 2014.
  • [24] NASA. "Six-Degree-of-Freedom Equations of Motion." Technical Reference.

Materials & Manufacturing

  • [25] Callister, William D. Materials Science and Engineering: An Introduction. Wiley, 2018.
  • [26] ASM International. "Aerospace Materials Handbook." 2021.
  • [27] Hexcel Corporation. "Carbon Fiber Composite Materials for Aerospace." Technical Data Sheet, 2023.
  • [28] Aircraft Spruce & Specialty Co. "Aircraft Materials Catalog." 2023.
  • [29] NASA. "Thermal Protection Systems for Hypersonic Vehicles." Technical Report, 2021.

Online Technical Resources

NASA Official Website
European Space Agency
SUPARCO (Pakistan)
Apogee Rockets Education
National Association of Rocketry
Aerotech Rocket Motors

Academic Disclaimer

This research project is conducted strictly for academic and educational purposes as a Final Year Project (FYP) at an HEC-recognized university in Pakistan. All information presented is based on publicly available sources and is intended for learning aerospace engineering principles. This project complies with PCAA and SUPARCO regulations for amateur rocketry and is NOT intended for any military or weaponization applications.

JOIN THE MISSION

Want to Contribute?

If you're an engineering student or enthusiast who wants to contribute to this project — purely for learning, gaining hands-on experience, and expanding your knowledge in rocketry, aerodynamics, or GNC systems — fill out the form below. This is a non-commercial academic project. I'll get back to you and we can discuss how you can be part of this mission.

Contact Form